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Neutral massive stable particles (dark matter candidates) are produced at colliders in pairs, due to 
an exact symmetry called a "parity" . They escape from detection, rendering their mass measurement 
difficult. We consider the pair production of these stable particles associated with two observable 
particles, from a two-step cascade decay of a heavier particle with even parity, via an intermediate 
particle with odd parity. We observe kinematic cusp structures in the invariant mass and angular 
distributions of the observable particles. Knowing the parent mass from direct resonant decay into 
Standard Model particles, one can determine the missing particle mass as well as the intermediate 
particle mass by using the cusped peak and end point of the distributions. The shape of the cusp 
distribution does not depend on the spin correlation. 

PACS numbers: ll.80.Cr, 12.60.-i, 14.80.-j 
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Introduction . Fault's postulation of a new particle that 
escapes from detection and carries away energy and an- 
gular momentum in f3 decay not only laid out the foun- 
dation for the weak interaction, but also rightfully in- 
troduced the first dark matter (DM) particle, the neu- 
trino. Ever since then, attempts to determine the masses 
and other properties of the neutrinos have led to many 
research efforts in nuclear physics, particle physics, as- 
troparticle physics and cosmology. If the upcoming ex- 
periments at the CERN LHC find evidence of large miss- 
ing energy events beyond the Standard Model (SM) ex- 
pectations, it will be a new exciting discovery. Most sig- 
nificantly, this may hold the key to explain the missing 
mass puzzle in the Universe, the dark matter. It is thus 
of fundamental importance to determine the mass and 
properties of this missing particle in LHC experiments, 
to uncover its underlying dynamics and to check its con- 
sistency with dark matter expectations. 

However, this task appears to be extremely difficult, 
even with the establishment of missing energy events at 
the LHC. In hadronic collisions, the undetermined longi- 
tudinal boost of the parton-level scattering leads to the 
ambiguity of their partonic cm. energy in the lab frame. 
Furthermore, with a conserved discrete quantum number 
(generically called a "parity" ) that keeps the lightest par- 
ticle in the new sector stable, the missing particles always 
come in pairs and thus the final state kinematics is even 
less constrained. Great efforts have been made to re- 
construct the mass of the missing particle in SUSY-like 
cascade decays. The mass difference of two parity-odd 
particles may be measured by the maximum end-points 
of invariant mass distributions [l| . Attempts have been 
made to determine the missing particle mass by the con- 
struction of new transverse variables [3, With long 
cascade decay chains, it may be possible to determine 
the masses through the on-shell mass constraints com- 
bining information from multiple events 

We propose a more favorable configuration in which 
the two missing parity-odd particles {X) come along with 
two visible SM particles (a), from the decay of a heav- 




FIG. 1: The "antler" decay diagram, a heavy particle (D) to 
two visible SM particles (a) and two missing particles (X), 
via two on-shell intermediate particles (B). 



ier parity-even particle {D) through intermediate parity- 
odd particles {B), as shown in Fig. [T] We dub this de- 
cay topology as the "antler" diagram. The advantage 
of antler decays is that once the parent mass (m^i) is 
known, the masses of the missing particle (mx) and the 
intermediate particle {rrig) can be determined by mea- 
suring the energy-momenta of the visible particles (a) 
without combinatoric complications. 

The antler topology is common in many scenarios with 
DM particle candidates. Familiar examples include the 
following theoretically well-motivated models: 



MSSM 0, 
Z' SUSY 0, 
UED 0, 
LHT i. 
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The precondition is that the mass of particle D is known 
a priori. This can be achieved since its even-parity al- 
lows its decay into two observable SM particles. In this 
regard; the antler topology is equally applicable to a lep- 
ton collider where the cm. energy is accurately known. 

Here, we explore two kinematic distributions: (i) Maa, 
the invariant mass of the two visible particles oi and 02, 
and {ii) cos Q, cosine of the angle between one of the two 
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Alaa [GcV] 

FIG. 2: Normalized difFerential decay rates versus the invari- 
ant mass Maa for various combinations of masses as given in 
Table ID The vertical lines indicate the positions of the cusps 
in each Maa distribution. 
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FIG. 3: Normalized differential decay rates versus cos O in 
the 73-rest frame (thin curves) and in the pp lab frame with 
y/s — 14 TeV (thick curves) . The parameters of MASS I and 
Mass III are given in Table HI 



visible particles (say oi) and the pair cm. moving direc- 
tion in their cm. frame. We show these distributions in 
Figs. [2] and [31 The distributions have, in addition to the 
end points, unique geometrical structures, the cusps (de- 
noted by vertical lines), which carry direct information 
on TUB and mx- 

A cusp is the point where two curves with different 
slopes meet. Cusps in the antler decay have outstanding 
merits for determining the missing particle mass: (?) a 
cusp is experimentally easy to identify due to its pointed 
feature; its sharpness depends on the mass spectra of the 
particles involved; (ii) looking for a cusp is statistically 
advantageous since it has large (in most cases, maximum) 
event rate; (in) cusps can determine both the masses 
of the intermediate particle B and the missing particle 
X; (iv) there is no combinatoric complication due to its 
simple decay topology; (v) spin correlations of the decay 
processes do not change the position of the cusps. 

The peak features of the cusps arise from the differ- 
ent origins for the Maa and cos0 distributions. Roughly 
speaking, the triangular shape of the Maa distribution 
comes from folding diagonally the flat distribution of 
d^r/d cos Videos 6*2, where 9i and 02 are the scattering 
angles of two visible particles. Reorganized along Maa, 
the two apexes of the triangle at both ends correspond to 
M™" and M^^, while the middle apex corresponds to 
the cusp. The largest overlapping area below the middle 
apex explains the largest phase space density. In the dis- 
tribution about cosO, the on-shell mass condition allows 
only a limited range for cos 8. Since the Jacobean factor 
leads to a steep increase at the end points, we have a 
cusp feature in the cos 8 distribution. We will leave the 
details to a future work 0. 

Cusp and Edge in Maa Distribution. We first only show 
the phase space distributions for on-shell particles. It is 
convenient to use the rapidities: the rapidity j] of particle 
B and the rapidity C of particle a in the rest frames of 
their parents D and B, respectively. In terms of the mass 



parameters, the rapidities rj and C are given by 
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Here and henceforth we use a shorthanded notation of 
Cx = cosh X. Obtaining the rapidities would be equivalent 
to measuring the masses ms and mx- 

(1) ma — case: Consider a to be massless first for 
simplicity. One would naively expect the invariant mass 
to have an end-point M^^ = mp — 2mx- However, due 
to the on-shell constraint for the particle B, we find a 
different end-point: 
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In addition, the Alaa distribution has a cusp at 



(2) 



(3) 



This is remarkable since the ratio MZ^^'/M^^^p = e^'' 
is governed by the initial decay D — > BB and thus gives 
m B ■ The product M™^^M^^^p depends on the secondary 
decay B — > aX and gives mx- Furthermore, dV/dMaa is 
of the form 



dV 



2r)Maa, 



if < Maa < M™«P; 
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Figure [2] shows dT/dMaa for four sets of representa- 
tive masses specified in Table ID The choice of the pa- 
rameters for Mass I is motivated by the Z*^^^ decay in 
the UED model Q • Since the two subsequent decays 

occur near the mass 
threshold, MASS I is to be called the "near threshold 
case" . For comparison, we consider the mass parame- 
ters with sizable gap in MASS III, the "large mass gap 
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mo (GeV) 


ms (GeV) 


ma (GeV) 


mx (GeV) 


Mass I 


1250 


600 





550 


Mass II 


1000 


440 





300 


Mass III 


1000 


350 





200 


Mass IV 


600 


250 


mz 


100 



TABLE I: Test mass spectrum sets for mass measurements 
using kinematic cusp structure (mz is the Z boson mass). 



case". The visibility of the cusp depends on the ratio 
Af^^'P/M^^''. As shown in Eq. (g]), the distribution for 
Maa < M^T^ is hnear, while that after M™''? is a con- 
cave curve with the maximum at M^'^/e. The cusp 
becomes a sharper peak if M^'^p > M^^/e (or equiva- 
lently tub > 0.44m£)). The parameters in Mass II are 
chosen to represent this boundary case of w OAAmu. 
The cusp structure is more pronounced for the near 
threshold case (MASS I), while it is less clear for the large 
mass gap case (Mass III). 

(2) iria 7^ case: If the SM particle a is massive (a Z 
boson or a top quark), we call it the "massive case" , as for 
the parameter choice in MASS IV, that are motivated by 
the MSSM heavy Higgs boson decay associated with two 
SM Z bosons. In this case, the analytic form of dT/dMaa 
is given by three pieces (the explicit forms are not very 
illuminating and thus not given here). The maximum of 

Maa is 



aa = 2maCn+C- 



(5) 



The positions of M™" and M^^'^p are as follows, depend- 
ing on the relations of the two rapidities rj and C,: 
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A few remarks are in order. First, the results of Eqs. ^ 
and ([6]) for 77 < C/2 include the massless case, which 
corresponds to C ^ oo- Second, for all three regions in 
Eq. ([S]), the Maa distribution shows a sharp cusp, as il- 
lustrated by the MASS IV case in Fig. [21 Third, the case 
C < ?7 has different Ai^™, rather than 2ma as naively 
expected. It is due to the enhanced boost of the two 
fast-moving parent i3's. This shift helps to resolve the 
ambiguity among the three regions. Still a two-fold am- 
biguity in the 77 < C/2 and (/2 < 7] < ( regions remains 
since we do not know whether the measured M^]^^p is 
2maCr,-Q or 2maCri- 

We propose another independent observable to break 
this ambiguity, (A|p!^|),„ax, which is the maximum of the 
difference between the magnitudes of the two transverse 
momenta of ai and 02. We find the expression to be 



max ( I ' 
ma [sinh(r/ ■ 



' C) - sinh \if] - C|] ■ 
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which is invariant under longitudinal boost. Note that 
the two-fold ambiguity happens when 77 < C leading to 
(A|p^|)niax = 2mo cosh C sinh 77. This sinh 77-dependence 
can provide independent information from M^^'p. 

Cusp in Angular Distribution. Another useful observable 
is the distribution about cos0, where Q is the angle of 
a visible particle, say ai, in the cm. frame of ai and 02, 
with respect to their cm. moving direction. We find the 
expression of dT /d cos O in the Z?-rest frame for rua ~ 
to be remarkably simple: 



dV 



dcosG 



sin O, if I cos8| < tanh77, 
0, otherwise. 



(8) 



The distribution has a sharp end-point, another cusp, 
with the highest event rate at the boundary: 



I cose I, = tanh77 = ^l-Am%/m 



(9) 



While this variable is unambiguous in the lab frame at a 
Icpton collider, wc cannot determine the longitudinal mo- 
tion of the particle D in hadronic collisions. After convo- 
luting with the parton distribution functions to convert 
it to the lab frame, dT/dcosQ in Eq. ([H]) is smeared. 
In Fig. [3l we compare the cos 9 distribution in the rest 
frame of D (thin curves) with that in the lab frame at the 
LHC (thick curves) for the near threshold case (MASS I) 
and the large mass gap case (MASS III). We have as- 
sumed that D is produced by direct s-channel gg or qq 
annihilation so that the longitudinal momentum distri- 
bution of D is obtained from the parton distribution of 
the incident protons. 

The convolution effects with the partons smear out the 
sharp cos cusp in the lab frame. For the near threshold 
case (Mass I), the two sharp rises at both ends get much 
less pronounced, although it is still possible to read the 
edge point off in the distribution. For the large mass gap 
case (Mass III), as the end point position approaches 
towards cos0 = ±1, the sharpness of two cusps main- 
tains better. It is interesting to note that the cusp in 
the Maa distribution and that in the cosO distribution 
provide complementary information for determining 771^ . 
The invariant mass distribution yields a better resolution 
for the near threshold case, while the angular distribution 
provides better one for the large mass gap case. 

To some extent, the cos 8 distribution in the D-rest 
frame may be obtained even in hadron collisions. The 
smearing effects can be effectively modeled by the well- 
known parton distribution functions in the large x region. 
Or the direct resonant decay of D into SM particles allows 
to extract the velocity distribution of D, which can be 
used to recover dF/dcosG in the D-rest frame. 

Discussion . Until now, we have ignored the effects of 
the matrix elements regarding the spin correlations be- 
tween the initial state and final state particles. We have 
confirmed that, for the four processes in Eq. ([T]), includ- 
ing the full matrix elements docs not change the shape of 
the distributions. The deviations from the phase space 
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Maa [GeV] 

FIG. 4: The invariant mass distribution with the effect of the 
finite decay width for the MASS I case. We vary the total 
decay width Tb with Tb = 0.01,5,60, and 300 GeV. The 
solid vertical line denotes the position of M^a^^. 

predictions become appreciable when the fermions and 
vector bosons (like the Z^^^^ in UED) have chiral cou- 
pHngs for both D BB and B aX decays. Even in 
this extreme case, the cuspcd peak remains at the same 
location and its height is changed by about 2%. 

However, non-vanishing decay widths of parent parti- 
cles can alter the shape of the distributions. In Fig. [U 
we show the Maa distributions for different finite de- 
cay widths Tb with the MASS I spectrum. With Tb ^ 
ms/lOO, such distortions do not have a big effect. With 
Tb ^ itlb/W, the finite decay width effect can change the 
cusp shape as well as the positions of M^^^p and M^^. 
In the limit Tb ~ tub , the momenta of the visible parti- 
cles oi and 02 span all of the allowed phase space given 
by mD and mx, regardless of tob, and thus the max- 
imum of Maa approaches mu — 2m x, denoted by the 
vertical dashed arrow in Fig. 5) The end point measure- 
ment with the known rriD leads to the missing particle 
mass mx^ just like the direct neutrino mass determina- 
tion in tritium decays. In the scenarios listed in Eq. ([T]), 



the intermediate particles (x2, L'^^^ and i_) typically 
have decay widths that are smaller than one percent of 
their masses. The cusp shape in the Maa distribution 
remains intact. 

Our cusp distributions of the signals arc quite distinc- 
tive from the SM background kinematics. For instance, 
the leading SM backgrounds for are mainly 

from W+W-/ZZ i+e'vD, while those for ZZ +$t 
from ZZZ/W^W~ Z , which should be well predicted in 
the SM. The optimal strategy to suppress the SM back- 
grounds is different case by case. We leave the detailed 
analyses to [1]. 

Conclusions . New techniques to measure the missing par- 
ticle mass are proposed based on an "antler decay topol- 
ogy" [D ^ B + B ^ aX + aX), with a final state of 
two visible particles (a) and two missing particles {X) 
via intermediate particles {B). This antler decay has no 
combinatoric complications. We found cusp structures 
in the invariant mass distribution of two visible parti- 
cles and an angular distribution. The cusp structure in 
the antler decay has large (in most cases largest) event 
rates unlike the end points. With its unique geometric 
feature, the cusp is easy to identify. Along with the end 
points of the distributions, the missing particle mass as 
well as the intermediate particle mass can be determined. 
A possible two-fold ambiguity can be resolved by another 
new observable, the difference between the magnitudes of 
transverse momenta of two observable SM particles. We 
believe that this technique will be invaluable for searches 
for new physics at the LHC. 
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